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A generalized analysis of the local entropy production of a simple fluid is used to show that, 
if intrinsic angular momentum is taken into account, rotational viscosity must arise in the linear 
non-equilibrium regime. As a consequence, the stress tensor of dense rotating matter, such as the 
one present in neutron stars, posseses a significant non-vansishing antisymmetrical part. A simple 
argument suggests that, due to the extreme magnetic fields present in neutron stars, the relaxation 
time associated to rotational viscosity is large (~ 10^^ s) . The formalism leads to generalized Navier- 
Stokes equations useful in neutron star physics which involve vorticity in the linear regime. 



1. INTRODUCTION 
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^ , In a recent paper |l[] , Rezania and Maartens have shown that relevant effects due to vorticity are present in neutron 
[t I ' star physics. Their communication was primarily based in results derived from kinetic theory which show that, 

^ , if Boltzmann equation is treated by means of Grad's method, a coupling between vorticity and shear necessarily 

' arises This paper is intended to present a study of an alternative role that vorticity may play in neutron stars. 
' This approach is based on a purely phenomenological formalism capable of describing the thermodynamical role of 
T-H ' vorticity and rotational viscosity, without introducing any assumption foreign to the concept of local equilibrium and 
^ ' the validity of the main conservation laws . 

A word of caution is here pertinent. In principle one could require the use of relativistic irreversible thermodynamics 
' to correctly cope with this question as it has been developed in previous work |3| . Nevertheless, in order to compare our 
: results with those obtained in Ref. ^ we shall overlook this fact and perform a non-relativistic analysis to emphasize 
on the role played by the rotational viscosity. 
7-H The paper is divided as follows: Section 2 gives a review of the Meixner-Prigogine approach to non-equilibrium 

thermodynamics Q. Attention should be paid to the fact that in this scheme, nothing forbids the existence of a 
non- vanishing antisymmetric part of the stress-energy tensor. Section 3 is devoted to the analysis of the constitutive 
^ ^. equations relating the thermodynamical forces and fluxes, and its immediate effects in neutron star dynamics. Section 
4 includes a brief comparison between the vorticity coupling suggested in Ref. and the one here presented. An 
outline of how it is possible to support the existence of a non-vanishing antisymmetric part of the stress tensor with 
kinetic theory is also included. 



2. MEIXNER-PRIGOGINE FORMALISM INCLUDING INTRINSIC ANGULAR MOMENTUM 
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j 2.1 Conservation laws 

The Meixner-Prigogine approach to irreversible thermodynamics is based in the consideration of the conservation 
laws of mass, linear momentum, angular momentum and total energy, plus the assumption of local thermodynamical 
equilibrium. The balance equations for mass and linear momentum, written in differential form, are respectively given 
by: 

f7 + (p-0, = §f + p(-0, = o (1) 

+ {pu'u'^ + r''^),^^ = p^+ rf/ = f (2) 
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where the summation is taken for repeated indices, p is the mass density, the local hydrodynamic velocity, r^'^ the 
stress tensor and the external force density term. The forces of interest in our formalism correspond to gravitational 
and magnetic fields. Latin indices will run from 1 to 3, and ; indicates a covariant derivative. 

Total angular momentum per unit of mass in any local element is properly described by an antisymmetric tensor 
jU _ rpj^-g ^^ej^gor possesses an orbital part L''^ and an intrinsic contribution 5*°^, so that 

jM ^ ^kl ^ gkl 

Since total momentum is conserved in purely rotating systems, we may write 

p^ + {x'r^'-x'r^').^^^0 (4) 

where x'' stands for the position vector. It is important to stress at this stage that the external force f'^ does not 
contribute to equation (^) for two reasons. The first obvious one is that the gravitational force being radial generates 
no torques. The second one is a little bit more subtle. Assuming that the external magnetic field is constant and 
points along the rotation axis, the Lorentz force generates only a radial component whence again generates no torque. 
On the other hand internal magnetic fields are accounted for in the intrinsic angular momentum equation as it will 
be pointed out below. Now, orbital angular momentum per unit of mass is given by 

L = X u — X u (5) 

From equations (|[) and (j^), it follows that, in the presence of the fields of interest, 

D Tkl 

SO that, substracting equation from equation we have, for the intrinsic angular momentum per unit of mass, 
that, 

p^^l^ = r'^^ - r" = -2r[«l (7) 
^ Dt ^ ' 

where square brackets on indices indicate the skew part. Intrinsic angular momentum, i.e. spin, in any local element 
of a neutron star is assumed to he proportional to its local magnetization M^^ , per unit of area, per unit of charge, so 
that: 

gki^pj^ke 

where the coefficient /3 has units of {Lenght)^ and is related to thermodynamical properties of the element. This 
coefficient will be estimated in Section 4. It follows then, from the corresponding balance of instrinsic rotational 
energy, per unit of mass, ej^ot], that using equations (0) and (^) 

0^ + {peirot]U = p-^ = — = -M,,rl (9) 

where the second equality is due to the fact that, the rotational contribution to the internal energy in a paramagnetic 
substance is proportional to B^^B^i and since the magnetization is proportional to the magnetic field S'"'', equation 
follows directly. 

Let us now consider the energy balance. Since total energy is conserved, one can write for this conservation law 
the expression, 



9 ( e[mec] + e[i 



dt 



{j[mec] + J[Q] + e[™t]w') ^ = (10) 



where Jjgj is the heat flux density and e[int\u^ the drift transport of internal energy. The explicit form of the heat 
flux JjQj will be given later. 

Multiplying equation (ph by the covariant vector ug and making use of equations (||) and (ph we get that: 



d e\ 



d 



f^ + (jU])^^ = ".;.-""M,,rM (11) 
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where 

1 



and 



e[mecl = 7:PU^ + Pf + P e[rot] (12) 



J[mec] = \pU^u'' + pipu^ + p e[rot]u'' + Ug t'^ (13) 



where ip stands for the gravitational potential and use has been made of the fact that no dissipation arises from the 
magnetic field, since ui e^.^ u^B^ = 0. If one now combines equations (|l|) and (|ll|), it is possible to obtain the balance 
equation for the internal energy, 



If the stress tensor is splitted into two parts, one related to the (scalar) equilibrium pressure p and one associated 
with the viscous pressure tensor H^*^ , then t^^ ^ p5^^ +11^'^, and equation ( p^ reads: 

D ^[int] 

Dt 

It must be noticed that rl^^'^l = nl'"''! 



{j[Q]) ^ - ut,k n"- - p «J + M.fcHl"! (15) 



3. LOCAL THERMODYNAMICAL EQUILIBRIUM AND ENTROPY BALANCE 

We now turn ourselves towards the question of introducing the entropy production. For this purpose we invoke the 
local equilibrium assumption which is crucial to the formulation of several versions of irreversible thermodynamics 
1^ . This assumption states that the local entropy s of the system is a time- independent functional of the local scalar 
thermodynamic densities, in this case the mass density p and the internal energy density e[j„t] (x*'',t). Then, 

s = s [p (x^ i) , e[™t] (x^^)] (16) 

whence, 

D s _ f d s\ Dp f ds \ De[int] 



(17) 



By the local equilibrium assumption per se, the thermodynamic coefhcients are the local expressions of their equilib- 
rium counterparts, 

^ P_ ( \ (m 

where p and T are the local pressure and temperature, respectively, and the time rates of change of p and e^nt] are 
given through equations (|I]) and ( p^ respectively. Putting all this information into equation (^^ and after some 
algebraic manipulation, one arrives at an expression of the form 

From standard irreversible thermodynamics, the entropy production or Clausius' uncompensated heat term E is 
precisely the right hand side of equation (|l9[) namely. 



^ = -^|§-^-.^n" + lM«n["i (20) 

Only when E > 0, which depends on the nature of the constitutive equations adopted, the theory will be consistent 
with the second law of thermodynamics. We will come back to this questions in the next section. 
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4. CONSTITUTIVE RELATIONS AND NEUTRON STAR DYNAMICS 



4.1 Linear constitutive relations 



We are now able to propose constitutive relations for thermodynamical fluxes and forces using Curie's principle [4, 
5]. This principle states that, given an expression such as (po| ) for the entropy production in an isotropic system, 
only tensors of equal rank may be coupled. According to this idea, we must decompose tensors H^'^ and ui-k in their 
irreducible forms, namely: 



and 



Ui-k ^ u[e-k] + U(e:k) + {u%) k 



(21) 



(22) 



where angular brackets denote the symmetric trace- free part and 11 represents one third of the trace of 11^'°. Intro- 
ducing these relations in equation (BO) we obtain: 



(23) 



Since we wish to have E > 0, linear constitutive relations are admissible to relate thermodynamical fluxes and forces, 
so that 



[Q] 



dT 

dxi 



(24) 



n 



(25) 



(l;k) 



and 



[l:k] 



M 



flk 



(26) 



(27) 



The transport coefficients introduced here are k, the heat conductivity, 77, the bulk viscosity, C, the shear viscosity, 
and rjn the rotational viscosity. The property sought for S will be satisfied if these coefficients are themselves positive. 
As is well known, this is borne out by experiment. 

4.2 Navier-Stokes equations in the presence of rotational viscosity and relaxation time. 



An immediate consequence of the validity of the constitutive relations ( |25[]27| ) is the effect of vorticity in the equation 
of motion (0), since 



d 

9^ 



n["i +n<"-> + (n + p) s- 



ek 



or 



Q^lk 

dx^ 



d 

dx^ 



dp 
dxk 



(28) 



(29) 



Now, as de Groot and Mazur observe this effect of vorticity involves a relaxation time obtainable from the 
assumption that ul^-'^l does not change substantially in a time interval < < < F where T is the relaxation time 
characterizing the local magnetization and defined below in equation (^3|). Thus, from equations (|27| ) and (Q), we 
have that, 



4 



2n[^ '=1 = (u^'-^ ''1 - M' = -P^^ (30) 

or, in terms of the local magnetization, equation (||) 

4r;H (.[^;'=l-M")=p/3:^ (31) 
which yields, under the uniformity assumption mentioned before, and assuming a vanishing magnetization at i = 0: 



(32) 



This coupling between intrinsic angular momentum and vorticity involves a relaxation time given by: 



r . f- (33) 

Densities in neutron stars are extremely large, typically near to 10^® kg m"^, so that T essentially depends on the ratio 
13/riR. Its magnitude will be explored below. 

4.3 Order of magnitude for the relaxation time 

The other parameter left in order to obtain a good idea about the relevance of this vorticity-magnetization coupling 
is fi/rjH in equation (|3^). Strong magnetic fields in the star, close to i? = 10^ Tesla, totally associated with local 
spin, suggest that the coefficient (3 may be estimated by a very simple argument, namely the relation between local 
magnetization, per unit of area, per unit of charge, and spin per kilogram of mass of a perfect paramagnetic solid, 
which is given by 

'"'ijcikm""' '''' 

where h is Planck's constant, /i = 10^^'* A.m^ is Bohr's magneton and C{T) is the Langevin function at temperature 
T. In our approximation, ^B/kT « 10^"''. Direct comparison with equation (^ yields, 

and, thus 

r.H^ (36) 

VR 

which may be quite significant if rja is a relatively small number. Estimates of this coefficient are scarce in the 
literature, but even if we assume that it is of the order of 10'^ Pa s^, F would be very large, if compared with 
relaxation times of the so-called "fast variables" associated with extended irreversible thermodynamics . 



5. FINAL REMARKS 



It must be stated that the formalism here developed is not in conflict with the one presented in Ref.[l], where 
vorticity appeared without taking into account intrinsic spin, but by means of a coupling between shear and vorticity. 
Whether the vorticity effects associated to intrinsic spin are relevant or not to stability of neutron stars is not clear. 
Boltzmann equation itself cannot predict values of rotational viscosity, since it leads to symmetric stress tensors and 
involves no information about internal degrees of freedom. A kinetic treatment involving the generalization suggested 
by Waldmann and Snider [9, 10] may be useful in order to compute precisely the value of rjn in equation ( ^3|) and so 
conclude about the real significance of this relaxation time. 

On the other hand, if the dynamical effect of spin were negligible in dense systems like neutron stars, then, according 
to equation (|3^), rotational viscosity must have large values and, from this point of view, further studies related to 
astrophysical implications of this coefficient should developed. 
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